Model study of adsorbed metallic quantum dots: Na on Cu(lll) 
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We model electronic properties of the second monolayer Na adatom islands (quantum dots) on 
the Cu(lll) surface covered homogeneously by the first Na monolayer. An axially-symmetric three- 
dimensional jellium model, taking into account the effects due to the first Na monolayer and the Cu 
substrate, has been developed. The electronic structure is solved within the local-density approx- 
imation of the density- functional theory using a real-space multigrid method. The model enables 
the study of systems consisting of thousands of Na-atoms. The results for the local density of states 
are compared with differential conductance (dl/dV) spectra and constant current topographs from 
Scanning Tunneling Microscopy. 



I. INTRODUCTION 

At certain faces of metals, such as the (111) face of 
noble metals, the surface electron states are confined to 
the vicinity of the top layer by the vacuum barrier x>n the 
vacuum side and the bandgap on the substrate sidetl. The 
electrons in these surfac^. slates form a two-dimensional 
nearly free electron gaaffl'tlB. It has also been observed 
that when adsorbing one to several monolayers of alkali 
atoms on these surfaces a manifold of discrete standing 
wave states, so called quantum welLstates (QWS), per- 
pendicular to the surface are formedBQ. These states can 
be detected, for instance, in photoemission spectroscopy 
(PES)@, inverse photoelectron spectroscopy (LEjES)H, two 
photon photoemission spectroscopyp(2PPES)E2land scan- 
ning tunneling microscopy (STM)til. A large amount 
of exper^gjjited j4ata is available for the system Na on 
Cu(lll)&QSBlia'y. The electronic structure and dynam- 
ics for this system have also beep-.kivestigated by first- 
principles theoretical calculationglSO. 

These localized surface states are of great interest since 
they play an important role in many ph^sieal processes 
like epitpdal growthtfJ, surface catalysigi3't3, molecular 
orderingiH and adsorptiontS. Experimental tools like 
STM and PES play an important role in the investiga- 
tions, since they enable spatial and spectroscopic resolu- 
tion of the electron states. 

One important discovery is the confinement of surface 
state electrons in so called quantum corrals. These man- 
made nanoscale structures are formed by deliberately 
assembling adatoms to enclosed structures by STME3. 
Due to the small size of the corrals, quantum effects are 
present, and both spatial and spectroscopic properties of 
the confined states can be studied experimentally. A nat- 
ural way of forming low-dimensional structures on metal 
surfaces is by controlled growth of epitaxial layers. With 
an appropriate choice of deposition and annealing tem- 
peratures small islands, so called quantjtHp dots (QD), 
with variable shapes and sizes may form£3. The advan- 
tage of these structures, in comparison with the corrals, is 
that they are relatively stable at low temperatures. This 



enables the imaging and investigation of their properties 
without inducing structural damage. One quantum me- 
chanical effect of the confinement is the increase of the 
surface state energy which in turn may lead to the de- 
population of the surface state band and thereby changes 
in the surface properties. 

In this paper we present calculated results for the elec- 
tronic structure of Na on Cu(lll), with the emphasis on 
describing the real-space resolved density of states nearby 
a sodium QD adsorbed on a sodium-covered Cu(lll) sur- 
face. Previously, an all-electron density-functional theory 
(DFT) study oLa free-standing Na layer in vacuum has 
been presentedEil. More recently, a DFT calculation us- 
ing ultrasoft pseudopotentials for the free-standing Na 
layer as well as-£pr the layers adsorbed on Cu(lll) have 
been presentec£3. A simple free electron model calcula- 
tion for a free-standing Na QD has_already been pub- 
lished by two of the present authorsE3. 

The calculations in this work are made in the context 
of the DFTBEi More specifically, the Rayleigh-qwjjtipnt 
multigrid (RQMG) method in axial symmetryOES is 
used for the numerical solution of the ensuing Kohn- 
Sham equations. The elect ratu ion interaction is simpli- 
fied using the jellium modelES, where the ions are re- 
placed by a rigid positive background charge of constant 
density. This model has provided basic physical un- 
derstanding of the electronic structures of simple metal 
surfaced23 thin metal filmsEj, vacancies and voids-inside 
metalsEa and finite clusters of simple metal atomsL2l. Re- 
cently, also uniform cylindrical aapowires have been stud- 
ied within the jellium modelBHLil 

The paper is organized as follows: Section [□] gives 
a short review over experimental results for the system 
Na on Cu(lll), Section III describes the computational 
method used in the calculations, Section |lV| discusses the 
details of our jellium model and in Section ^ the results 
are presented and comparisons are made with experimen- 
tal findings. Finally, Section VI gives the conclusions. 
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II. NA ON CU(lll) 

Alkali metals adsorbed on the closed-packed (111) sur- 
face of metals form hexagonal structures at saturated 
monolayer coverages, following approximately the under- 
lying substrate structured. The first monolayer of Na 
on Cu(lll) is observed to saturate at the coverage of 
6 = 4/9 ~ 0.44E5I, corresponding to 4 Na atoms per 9 
surface Cu atoms. The Na atoms thus form a hexagonal 
(3/2 x 3/2) structure and the Na atom spacing of 7.43 ao 
is comparable to the atomic distance of 6.92 ao in bulk 
Na. 

The adsorption of Na atoms on the Cu(lll) surface 
will induce a charge redistribution at the interface be- 
tween the adlayer and the substrate. It has been seen 
from photoemission experimentst3&£l that when the Na 
coverage is increased the Cu Shockley surface state de- 
creases in energy. For coverages above O w 0.11 the sur- 
face state is shifted below the lower band edge of the local 
band gap of the Cu(lll) surface, and is no longer visible 
in photoemission-experiments. Two-photon photoemis- 
sion experimentslHj indicate other unoccupied Na-induced 
states in the local bandgap at the Cu(lll) surface, which 
will also decrease in energy with increasing Na coverage. 
For higher coverages, the lowest of these states will be 
down-shifted below the Fermi energy, and thus get occu- 
pied. At the saturated rpepejlayer coverage, this state will 
be located about 0.1 eVli2Hi3t^ below the Fermi energy at 
the T-point of the surface Brillouin zone (SBZ). The cor- 
responding next-lowest state will be located 2.1 eV above 
the Fermi levellj at the T-point. The lower of the two 
states has one node in the z-direction, while the higher 
has two. 

If the Na atom deposition continues after the first 
monolayer (ML) is completed, a sepapd layer will start to 
grow. Recent STM measurement sOH indicate that the 
second monolayer of Na grows via the formation of com- 
pact islands with hexagonal atomk arrangement. Nor- 
mal photoemission experimentsBEl indicate that when 
the second monolayer grows the emission intensity due 
to one-monolayer states decreases gradually, and for cov- 
erages above 1.3 ML a new peak, approximately 0.1 eV 
above the Fermi energy appears. This peak is ascribed to 
the two- monolayer thick parts, and the energy is shifted 
to somewhat lower values as the coverage is increased. 

Some theoretical attention has also been paid to Na 
on Cu(lll), including the island growth. Free-electron 



model calculations have bepj^erformed ictr circular NaE2l 
as well as hexagonal Agocil and Nac3 free-standing 
islands. All-electron calculations for an unsupported 
monolayer of NaEH and first-principles slab calculations 
for one atomic Na layer in (12.x 2) and (3/2 x 3/2) adsor- 
bate structures on Cu(lll)E3 have been presented. In the 
present paper, we report jellium model calculations for 
an unsupported monolayer of Na and a cylinder shaped 
free-standing Na QD. We also present two- density- jellium 
calculations for the system Na on Cu(lll), where we 
have modeled the underlying Cu(lll) substrate by us- 
ing a lower density slab to mimic the decay of the sur- 
face states into the substrate. Comparison is made with 
experiments and previous theoretical calculations. 



III. COMPUTATIONAL METHODS 

In the Kohn-Sham scheme of DFT, one solves the elec- 
tron density n(r) of the system selfconsistently from a set 
of equations. One of these equations is the single parti- 
cle Schrodinger equation. The models used in this work 
are axially symmetric. Thus the Schrodinger equation 
is separable, and the wave functions can be written as 
products 



tkn(r, z, 4>) = e im0 Z7 mkn (r, z) 



(1) 



Above, m is the azimuthal quantum number implied by 
the axial symmetry while n differentiates between orthog- 
onal states with same m and k. In the calculations in- 
volving the infinite monolayer, two k- vectors are used as 
explained in Sec. |y|. The external potential of the sys- 
tems studied in this work is caused by the positive back- 
ground charge n + (r) . The effective potential V e s includes 
also the Hartree potential of the electron density and the 
exchange-correlation potential Vxc ,r-which we treat in the 
local density approximation (LDA)c2l. The electron den- 
sity n (r) is obtained by summing single-electron densities 
with the occupation numbers / m kn- The degeneracies of 
the states are taken into account by the factor (2 — 5o m ) 
and the occupation numbers fmkn obey the Fermi-Dirac 
statistics with a Fermi level (Ep) so that the system is 
neutral. A finite temperature of 1200 K is used to sta- 
bilize the solution of the set of equations. Thus, in the 
present axial symmetry (r = (r, z)) the Kohn-Sham equa- 
tions read as 



1/19 d 2 m 2 d 2 



2 V r d~r + W 2 ~ ^d^ 2 + 2Veff ) = e ^U m i^ n (r) (2) 

n(r) = 2 J2( 2 ~ ^0m)/mkn|/7 mkn (r)| 2 , (3) 

mkn 



V eS (r) = <S>(r) + V xc (r), 



(4) 
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The Schrodinger equation is solved using the 
Rayleigh-quotient multigrid (RQMG) methodE3 which 
has been implemented— In various geometries, including 
the axial symmetryE§E3. In the RQMG method, the 
Rayleigh quotient (i/j\H\iJj) / (V^) on the finest level grid 
is directly minimized, the orthogonality constraint being 
taken into account by a penalty functional. The Poisson 
equation (||) is solved for the electrostatic potential 3>(r) 
using a standard multigrid methooO. 

To obtain selfconsistency, we use the simplest possible 
potential mixing scheme, 



(l - A)V? r 



(6) 



The largest system of this work contains 2550 electrons, 
the diameter of the supercell being 170 A. Obtaining self- 
consistency in such a system requires a very small A value 
of 0.005. Otherwise the charge sloshing results in diver- 
gence. More sophisticated mixing strategiesEZI will be in- 
dispensable in the future calculations. However, because 
of the simplicity of our model systems, we can very ac- 
curately estimate an initial guess for the selfconsistent 
effective potential of large systems using the more easily 
convergent smaller systems as reference. 

In our largest calculation a grid of 319 x 95 points is 
used for the presentation of the wave functions, potential 
and density. Taking into account the unoccupied states 
needed in the modeling, up to 2400 different states have 
to be solved at every selfconsistency iteration. Luckily, it 
is straightforward to parallelize the calculation over the 
65 different m values and the.-fcwo k-points (see below). 
Moreover, the RQMG methodE3 handles this part of the 
calculation with optimal efficiency. 

The two-jellium model for the surface described in Sec. 
IV results in an asymmetric density distribution with a 
surface dipole. Thus the electrostatic potential on the 
substrate side is a constant different from that on the 
vacuum side. For the Poisson equation, we thus use the 
boundary condition of zero derivative on the substrate 
side, and that of zero value on the vacuum side. Solving 
the Poisson equation, the boundaries above and below 
the system are extended to a distance five times greater 
than in the case of the wave functions. 

In this work, we calculate the local density of states 
(LDOS) above a surface at distances corresponding to 
those typical in STM measurements (of the order of 
20 ao). At such distances, the amplitude of the wave 
function decreases by several orders of magnitude. This 
kind of modeling is thus a serious test for the RQMG- 
method. We have checked the accuracy of our method for 
the spherical harmonic oscillator and a model hydrogen 
atom potential, for which the wave-functions are known 



analytically. The evanescent tails of the wave functions 
solved with the RQMG method agree with the analyt- 
ical ones, even when the amplitude of the wave func- 
tions has dropped by 20 orders of magnitude. This level 
of accuracy is beyond the reach of plane- wave methods, 
where periodic boundary conditions are necessary, and 
which provide a uniform accuracy across the calculation 
volume, resulting in spurious oscillations in the vacuum 
parts of the system. 



IV. MODELING THE SYSTEM 

We are interested in the system of a monolayer-thick 
Na QD on the complete Na monolayer on the Cu(lll) 
surface. We know from experiments that these islands 
are approximately hexagonal in shape, following the un- 
derlying structure of the Na monolayer. 

In order to interpret recent STM data for these types 
of systems, mapping the energy resolved real- space elec- 
tron density near a QD is necessary. First, it is of in- 
terest to find what level of theoreticaLjpadeling is re- 
quired. It has been shown previouslyCj'tUcS that sim- 
ple two-dimensional 'particle-in-a-box' calculations give 
qualitatively good results, in the sense that the peak 
structure of LDOS resembles spectra obtained in the 
STM dl/dV measurements. In this work we improve 
the theoretical description by performing selfconsistent 
three-dimensional DFT calculations, where the effects of 
the underlying monolayer and substrate are introduced. 
The hexagonal QD is modeled by a cylindrical jellium 
QD, and the underlying Na monolayer and Cu(lll) sub- 
strate by the two-density jellium slab, as described below. 
Comparisons between calculations of free-standing QD's 
and QD's on a substrate show indeed that the underlying 
monolayer and substrate induce a new type of states that 
the simple 'particle-in-a-box' calculations cannot account 
for. Since the z-dependence of wave- functions is included 
in our calculations, we can calculate the tunneling cur- 
rent at realistic STM-tip distances above the system, and 
estimate the energy dependence of the step height from 
the calculated constant current topographs. 

The different model systems studied in this work are 
shown in Fig. [l]. Our model is readily applicable to the 
case of a free-standing cylindrical quantum dot, where 
we use zero Dirichlet boundary conditions for the wave 
functions and for the Coulomb potential. The next step 
is to model a free-standing monolayer. A uniform planar 
system cannot be exactly reproduced in the axial sym- 
metry. We adopt an appmximation scheme analogous to 
the Wigner-Seitz methods. We imagine the plane being 



4 



a) 



b) 



c) 



d) 



FIG. 1: Profile of the axially symmetric background charge 
in the case of a) jellium model for free-standing Na quantum 
dot b) jellium model for free-standing Na monolayer c) two- 
jellium model for Na monolayer on Cu(lll) d) two-jellium 
model for Na quantum dot on Na monolayer on Cu(lll) 



filled by hexagons, and then approximate these hexagons 
by area-covering circles. In order to sample the Bril- 
louin zone of the lattice of circles we use two k-points, 
k = and k at the Brillouin zone boundary. The wave 
functions with k = are required to have a vanishing 
radial derivative at the radius of the circle whereas the 
wave functions with k at the Brillouin zone boundary 
vanish there. According to our calculations the model 
gives a uniform (r-independent) charge distribution for 
the monolayer. It also minimizes the interactions.. be- 
tween a QD inside a circle with its periodic imagescj. 

The next step is to place the Na monolayer on top 
of the Cu(lll) substrate. The effect of the substrate is 
modeled using the two-jellium model which is illustrated 
in Fig. |l|c and Fig. ||. We do not model the electrons 
of the bulk Cu. The density of electrons per unit area 
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FIG. 2: One complete monolayer of Na on Cu(lll) sur- 
face within the two-jellium model. The positive background 
charge (shaded areas), electron density (dashed line), effective 
potential (solid line), and electrostatic potential (dash-dotted 
line) are shown. The shading corresponds to Fig. |l|. 



in the two-jellium model is kept the same as in the jel- 
lium model for a free-standing monolayer. We add a 
layer of lower density jellium, in order to mimic the dif- 
ferent wave function decays into the substrate and into 
the vacuum. The thickness W2 and density (via r S 2) give 
two free parameters of the lower density jellium, which 
we adjust in order to reproduce the relevant experimen- 
tal values of the first and second surface band bottoms 
at the coverages of 1 ML and 2 ML, respectively. Here 
the first and second bands correspond to wave functions 
with one and two nodes in the vertical direction, respec- 
tively. The thickness w[ of the higher density jellium in 
the two-jellium model is given by 

w[ = Wi - ( — ) w 2 , (7) 

where w\ is the thickness of the free-standing Na mono- 
layer. 



A. The underlying monolayer and substrate 

To test our model, we first study the systems of a free- 
standing Na monolayer and that of a Na monolayer on 
Cu(lll), and compare the results with other theoretical 
results and experimental findings. The Na jellium density 
is determined from the bulk nearest neighbor distance of 
6.92 ao and the experimenjk&L-height of 5.5 ao (2.9 A) 
of 1 ML of Na on Cu(lll)laM The resulting density 
parameter, r s = 3.79 ao, gives a slightly higher density 
than its bulk value of 3.93 ao for Na. The thickness and 
the density of the lower density slab have been chosen by 
fitting the bottom of the second band for the 1 ML Na 
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FIG. 3: The uppermost panel: Energy spectrum of the free- 
standing cylindrical jellium QD (Fig. fila) containing 100 elec- 
trons. The discrete eigenenergies are shown as a function of 
the quantum number m. The top of the figure corresponds 
to the vacuum level. The middle panel: LDOS calculated at 
the cylinder axis at 8 ao above the jellium edge. The lowest 
panel: LDOS calculated at the cylinder axis at the jellium 
edge. 



coverage and that of the third band for the 2J dL . |Na cov- 
erage on Cu(lll) to the experimental valuei3c3Ej. The 
values of r S 2 = 6.0 ao and W2 — 6.3 ao give (using the 
unit cell of radius 72.8 ao containing 400 electrons per 
monolayer) in the 1 ML case the bottom of the second 
band at 75 meV below the Fermi level and in the 2 ML 
case the bottom of the third band at 50 meV above the 
Fermi level. These values are reasonably close to the ex- 
perimental values of abojot-100 meV below and above the 
Fermi level, respectivelyBl3L2l. The correct positions rel- 
ative to the Fermi level are important, because we solve 
for the electronic structures selfconsistently, so that the 
occupancies of the single-electron states affect the poten- 
tial and the character of the states themselves. 
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FIG. 4: The uppermost panel: Energy spectrum of the sys- 
tem of a QD containing 100 electrons on top of a two-density 
jellium slab described by a supercell of 400 electrons (Fig. |l|d). 
The k = eigenenergies are given by by thick horizontal bars. 
The thin vertical bars indicate the dispersion in the k- space. 
The middle panel: LDOS calculated at the cylinder axis at 
14 ao above the jellium edge. The lowest panel: LDOS calcu- 
lated at the cylinder axis at the jellium edge. 



B. The quantum dot 

We start by studying a free-standing monolayer-thick 
Na QD (Fig. 0a), since this system shows close resem- 
blance to the simple particle-in-a-box system often used 
as a first approximation when describing the electronic 
structure of a QD on a surface. The number of atoms 
in the QD is chosen to 100, which corresponds to a QD 
radius of about R = 36.40 ao- The uppermost panel in 
Fig. H shows the corresponding energy spectrum relative 
to the Fermi energy Ep. Note that the the discrete en- 
ergy eigenvalues are plotted as a function of the quantum 
number m and not as a function of k. There are three 
bands below the vacuum level, but only the first band 
(no horizontal nodal planes) is occupied. The emergence 
of the succeeding second and third bands can be seen as 
the condensation of the energy levels at around 1.1 eV 
and 2.5 eV, respectively. 

The LDOS calculated at the cylinder axis at the jel- 
lium edge and at 8 ao above the edge are shown in the 
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FIG. 5: Wave functions for the 21 lowest lying m — 0, k = 
states for the system of a monolayer-thick QD, containing 100 
electrons, on top of a Na/Cu slab (400 electrons per super- 
cell) . The wave functions are plotted in a plane parallel to the 
z-axis through the center of the QD. In each subfigure, the 
cylinder axis is shown by the solid vertical line. The shading 
indicates the positive background charge and the dashed ver- 
tical line points its QD edge. The upper, lower and right-hand 
subngureborders and the cylinder axis limit the computation 
volume with the dimension of 73x60 ao 2 . 



lowest and the middle panel of Fig. 3, respectively. Only 
states with m = contribute, since they are the only ones 
with nonzero contributions at the axis. The discrete en- 
ergy levels are broadened to Lorenzians with the width 
T = 8 meV. The LDOS at the jellium edge can easily 
be resolved in terms of the contributions from the differ- 
ent bands: The peaks corresponding to first, second, and 
third bands form series with quadratically increasing in- 
tervals and smoothly increasing peak amplitudes. At the 
distance of 8 ao above the QD, the contribution due to 
the first band states is diminished and the contribution 
due to the third band states with high quantization in 
the z-direction (two horizontal nodal planes) is dominat- 
ing the LDOS. Comparison with the results of the simpjkj 
'particle-in-a-box' calculation by Lindberg and HellsingES 



shows that these two calculations give qualitatively the 
same results. 

We now compare these results for the free-standing 
QD to the system of the QD adsorbed on a Na mono- 
layer on the Cu(lll) surface (Fig. |l|d). In our calcula- 
tion, the substrate is a two-jellium cylindrical supercell 
containing 400 electrons. The energy spectrum is shown 
in the uppermost panel of Fig. |] in which the different 
m levels corresponding to the k = -points are given 
with the k dispersion calculated using the k point at the 
Brillouin zone boundary. In comparison with the spec- 
trum of the monolayer-thick QD in Fig. || the bands are 
shifted downwards, because of the larger jellium thick- 
ness at the QD. The lowest energy states (in the bulb 
of the level diagram) have no k-dispersion and they are 
localized at the QD and the substrate slab below it (See 
Fig. |^). Their dispersion as a function of m is similar 
to that in Fig. [|. Introducing the underlying substrate 
gives rise to a new type of states, which are not localized 
to the QD region. In fact, these states form the over- 
whelming majority. As a result, new bands are induced 
in the energy spectrum and they are less dispersive as 
a function of m. The reduced m dispersion reflects the 
fact that the states are extended over the entire circular 
supercell. For the same reason, these bands have a larger 
dispersion in the k-space than the localized QD bands. 
A simple particle-in-a-box or free-standing QD calcula- 
tion cannot provide states of this kind, and it is therefore 
interesting to see to what extent they contribute to the 
local electronic structure above the QD. 

The LDOS for the QD adsorbed on a Na monolayer on 
the Cu(lll) surface is given in the lowest and the middle 
panel of Fig. ^ at the jellium edge and at 14 ao above the 
edge, respectively. In order to avoid complications due to 
the interactions between the supercells we calculate the 
LDOS using only the k = states and the LDOS is then 
calculated as in the case of the free-standing monolayer- 
thick Na QD. To enable a thorough comparison with the 
free-standing QD results we have to study first the wave 
functions in more details. 

Fig. U gives all the wave functions in the interesting en- 
ergy region for the m = states of the QD adsorbed on 
the Na monolayer. The first three states, n = 1, 2, and 
3, correspond to states within the first band. They are 
localized to the QD and the substrate slab below it and 
they have no nodes in the z-direction. The n = 4 state 
shows another character with a density no longer local- 
ized to the QD region but spread also to the slab region 
around. This is the first state belonging to the new type 
of bands, induced by the slab. State n = 8 is the begin- 
ning of the next band consisting of states with one node 
in the z-direction (second band in the QD). These states 
are resonance states, the amplitude of which is strongly 
enhanced in the QD region, but due to the hybridization 
with the delocalized slab states they are actually delo- 
calized to the whole system. The n = 14 state starts the 
next band consisting of delocalized states with one hori- 
zontal node in the slab region, i.e. it is a second band in 
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the slab. Finally, state n = 16 represents the first reso- 
nance state with two nodes in the z-direction (third band 
in the QD). In Fig. || one notes that the states n = 16 
and 17 and also the states n = 19 and 20 form pairs. The 
state lower in energy in the pair is a bonding combination 
of a QD state and a surrounding slab state whereas the 
state higher in energy is an antibonding combination. 

The electronic structure of the QD adsorbed on the Na 
monolayer, discussed above in terms of the m = wave 
functions, is reflected in the LDOS in Fig. [|. First, in the 
LDOS at the jellium edge we notice that, in comparison 
with the free-standing QD model, the underlying slab 
both introduces new type of bands and squeezes the other 
bands more tightly together to fit more states below the 
Fermi level (the lowest panels of Figs. 3 and 4). There- 
fore the present LDOS looks qualitatively different from 
that for the monolayer-thick unsupported QD in Fig. [|. 
Moreover, the hybridization of the QD and surround- 
ing slab states to bonding-antibonding pairs causes the 
splitting of peaks seen clearly for the third band states at 
the distance of 14 ao above the jellium edge (the middle 
panel of Fig. 4). If we had a continuous spectrum of 
slab states we would have a single resonance peak with a 
finite energy width. 

Having the modeling of the STM results in mind, the 
interesting question arising is whether or not the localized 
states calculated by the free-standing QD model give the 
same LDOS far above the QD as the states in the model 
including the substrate slab. Studying the LDOS plots, 
matching each peak with the corresponding wave func- 
tion, we notice that the resonance states with strongly en- 
hanced amplitude in the QD region are dominant at large 
distances above the QD. The contribution of the more de- 
localized slab states is small. Therefore the free-standing 
QD model is expected to preserve validity in predicting 
LDOS at large distances above the QD. The too broad 
energy spectrum in the free-standing QD model can be 
corrected for by increasing the dot height with a mono- 
layer of Na jellium or with a two-jellium layer. 



V. COMPARISON WITH EXPERIMENT 

One of the most useful instruments in surface, science 
is the Scanning Tunneling Microscope (STM)c3'lZI. It can 
be used for measuring the real-space electronic distribu- 
tion with atomic resolution, as well as the local energy 
distribution of electrons and the lifetimes of excited elec- 
tron states. The real-space distribution is achieved when 
scanning the surface either in the constant current mode, 
where the tunneling current is kept constant by changing 
the tip-surface distance using a feedback mechanism, or 
in the constant height mode, where the tunneling current 
is measured when scanning the surface at a constant tip 
height. The resulting image then displays the topogra- 
phy of the surface. Information about the local electronic 
structure is obtained by measuring the current variation 
with the applied voltage. This quantity, the differential 
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FIG. 6: Cylindrical QD containing 550 electrons on two- 
jellium substrate. The local density of states is shown at 18 ao 
above jellium edge at the axis (solid line) and at r = 20 ao 
(dashed line) away from the axis. Lorenz broadening with 
r = 8 mpV has been used. The relative experimental peak 
posit ions£3 are given by vertical arrows pointing downwards. 
The peaks are identified with (m, N) resonance states having 
two horizontal node planes in the QD. 



conductance dl/dV, is proportional to the product of the 
local densitYj-of states (LDOS) and the transmission co- 
efficient XoEj. However, if the applied voltage is small, 
the bias dependence of T is small, and with Eq. || we 
approximately have 



dl 
dV 



oc 



^(2-(5 0m )|/7 mkn (r)| 2 ^ 



mkn 



iV). (8) 



mkn 



In the STM study by Kliewer and BerndtEl, constant cur- 
rent topographs and dl/dV measurements are presented 
for a Na island on Na monolayer on Cu(lll). The size 
of the island is 230 x 170 a 2 (120 x 90 A 2 ). We have 
studied a cylindrical jellium dot with similar dimensions, 
i. e. having the radius of 85 ao and containing thus about 
550 electrons. The Na/Cu substrate is described in our 
calculations by a cylindrical two-density-jellium supercell 
with the radius of 160 ao and containing 2000 electrons. 
The radius of 85 ao is actually fixed to reproduce the .peak 
structure of dl/dV spectra by Kliewer and Berndtc3 as 
well as possible (See Fig. || and [discussion below). For 
comparison, Kliewer and Berndtc3 used in their model- 
ing two-dimensional hard-wall hexagons with the radius 
of 83 ao- 

We show in Fig. || the LDOS at 18 ao above the jel- 
lium edge both at the z-axis of the QD and at r = 20 ao 
away from the axis. The height corresponds to a typ- 
ical tip-sample distance in the STM experiments. The 
LDO S is calculated as for the smaller systems in Sec. 
IV B . The peaks in the figure correspond to states with 
two horizontal nodes in the 2 ML part of the system 
(third band in the QD). At the axis, only the m = 
states contribute, while away from the axis also peaks 
with m/0 occur. The LDOS peaks can be labeled with 
the "quantum number" N by counting for the number 
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FIG. 8: Cylindrical QD containing 550 electrons on two- 
jellium substrate. The step height of the second Na mono- 
layer determined from calculated constant current surfaces 
(see Eq. ^ is shown as a function of the bias voltage (energy 
relative to the Fermi level). The inset shows the LDOS isosur- 
face profiles (height-to-radius ratio exaggerated) at energies 
-400 meV (dashed line) meV (solid line) and 400 meV 
(dash-dotted line). The height is measured from the jellium 
edge of the second ML QD. 
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FIG. 7: Cylindrical QD of 550 electrons on two-jellium sub- 
strate. Isosurfaces of the electron density (top) and the LDOS 
(with r = 0.8 meV) at energies corresponding to the domi- 
nant peaks of Fig. o. The quantum numbers of the dominant 
states contributing at each energy is indicated. The isovalue 
for each plot is chosen as the value of the corresponding quan- 
tity at 18 ao above the jellium edge and30 ao off from the axis. 
The height-to-radius ratio in the plots is exaggerated. 



of radial nodes of the corresponding wave functions in 
the 2 ML part (See Fig. |). For ra=0, the N = 1 state 
has no radial nodes in the dot region whereas N = 2 
has one radial node and so on for larger N. As in the 
case of the smaller dot discussed above in Sec. IV B, the 



states strongly peaked in the QD are resonance states due 
to the hybridization with the states of the surrounding 
monolayer and span the whole system. Besides the de- 
localization of the states, the resonance character causes 
the fact that in the LDOS (Fig. |6|) several peaks may 
correspond to the same resonance state as discussed ear- 
lier for the smaller system. We have identified the LDOS 



peaks by examining the wave functions. The horizontal 
lines below the quantum numbers m and N connect the 
peaks belonging to the resonance in question. 

The relative positions of peaks appearing in the ex- 
perimental dl/dV spectra by Kliewer and Berndtc3 are 
shown in Fig. || as arrows pointing downwards. The 
experimental data is shifted so that the lowest experi- 
mental peak coincides with the lowest calculated peak. 
The experimental spectrum is recorded slightly off from 
the center of the hexagonal QD which should be taken 
into account when comparing with the calculated results. 
Previous two-dimensional free-electron calculations for 
hexagonal potential boxes-Jiaye reproduced well the ex- 
perimental peak positionscil'El. In the present modeling, 
the experimental peak positions agree with the calcu- 
lated m = resonance positions with the exception of 
the third and fifth experimental peak. The calculated 
m = resonances obey the pattern s n = Eq + AN 2 , 
as would be expected for a free particle in a hard-wall 
cylinder. We have fitted E « 22 meV, A « 15 meV. It 
is gratifying to note that in the LDOS recorded off the 
cylinder axis strong m = 1 resonance peaks appear so 
that the third and fifth experimental peaks can be ex- 
plained. Thus, our model can reproduce quantitatively 
the experimental peak positions. According to our cal- 
culations the resonance width increases toward to higher 
energies. The increase is maybe slightly stronger than in 
experiment, indicating a somewhat too weak confinement 
of the resonance states in our model. 

We have also calculated the isosurfaces of the LDOS 
at the energies corresponding to the dominant peaks in 
Fig. |[ The results are shown along with the total elec- 
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tron density in Fig. [7[ In order to see clearly the nodal 
structures of the different states, the LDOS is calculated 
using a smaller lorenzian width of 0.8 meV. The density 
is smooth in the interior of the QD and shows minor os- 
cillations at the perimeter of the QD. The development 
of the nodal structure is clear and compares qualitatively 
well with that found in the experimental dl/dV mapsESl. 
It can also be seen that the isosurfaces corresponding 
to the two (m,N) = (0,4) peaks differ from each other 
mainly near the perimeter of the dot. The higher peak 
shown at the energy corresponding to the highest (0,3) 
state, on the other hand, does not show equally clear 
(0,3)-character. The explanation is a state with quan- 
tum numbers (2,2) at almost exactly the same energy. 

From Eq. || we obtain a simple formula for the tunnel- 
ing current, 

rE F + U 

I(U, r,z)oc / LDOS(£, r, z)dE. (9) 

JEf 

From these numerical constant current topographs, we 
can estimate the apparent step height at the perimeter 
of the QD. The inset of Fig. || shows the constant LDOS 
height as a function of the distance from the cylinder axis 
for the bias voltages of -400, 0, and +400 meV. The value 
of the LDOS is 10 -11 arb. units on the scale of Fig. ||. 
The absolute height of the isosurface from the jellium 
edge depends naturally on the LDOS value chosen, but, 
according to our calculations, the relative changes are in- 
sensitive to the LDOS value over a wide range of values. 
In order to construct the apparent step height we first 
obtain the numerical constant current topographs (Eq. 
|9|), then average the profiles over the oscillations above 
the 2 ML and 1 ML parts of the system and take the 
difference. The results are shown in Fig. || as a function 
of the bias voltage. The trends seen can be explained by 
studying the LDOS-isosurfaces (see the inset of Fig. ||), 
and then noting that the constant current surfaces are 
obtained by simple integation (Eq. ^). At -100 meV, the 
second band starts to contribute in the 1 ML part of the 
system rising the height there and thereby lowering the 
step height. Then the onset of the third band in the 2 
ML part rises the step again. This rise is similar to that 
seen in the experiment by Kliewer and BerndtH as well 
as the decline at higher voltages. However, the compar- 
ison with the experimental result shows differences: our 
step height is too low by a factor of 2, and the raising 
of the step at negative bias voltages is not seen in the 
experiment. There may be several reasons for the differ- 
ences in the step heights. One is that the experimental 
step height of 5.5 ao (2.9 A), which is determined at a 



voltage just before the rise in the step height, is directly 
used as the thickness of the jellium describing the second 
monolayer of Na. A more consistent procedure might be 
to take the voltage dependence into account. Moreover, 
the apparent step height depends on the relative vacuum 
decay rates of the second and third band states of the 1 
ML and 2 ML systems, respectively. Their correct de- 
scription may be too demanding for our simple model. 
VI. CONCLUSIONS 

In this paper we have presented a model for the elec- 
tronic structures of alkali metal islands or quantum dots 
adsorbed on metal surfaces. In particular, we have fo- 
cused on the system of Na on the Cu(lll) surface, where 
approximately hexagonal Na quantum dots have been ob- 
served to form during the epitaxial growth of the second 
Na monolayer. 

We have modeled the quantum dots as small cylindri- 
cal jellium islands, and the underlying Na monolayer and 
Cu substrate as a two-density jellium slab. The param- 
eters of the model have been chosen to fit experimental 
spectroscopic data and calculated first-principles band 
structures for one and two completed monolayers of Na 
on the Cu(lll) surface. The calculations were performed 
in the context of the density-functional theory, using a 
real-space electronic structure calculation method. 

The calculated results are compared with experimental 
findings from scanning tunneling microscope and photoe- 
mission experiments. The model gives local densities of 
states which are in a quantitative agreement with con- 
stant current topographs and dl/dV spectra and maps. 
Thereby the idea of surface states which are localized as 
resonances at the quantum dots is supported. The fu- 
ture applications of the model will include studies of the 
adsorption and dissociation of molecules in the vicinity 
of alkali metal quantum dots. 
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